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ABSTRACT 

Helium that accretes onto a Carbon/Oxygen white dwarf in the double white dwarf AM Canum Venaticorum 
(AM CVn) binaries undergoes unstable thermonuclear flashes when the orbital period is in the 3.5-25 minute 
range. At the shortest orbital periods (and highest accretion rates, M > 10~ 7 M Q yr" 1 ), the flashes are weak and 
likely lead to the Helium equivalent of classical nova outbursts. However, as the orbit widens and M drops, 
the mass required for the unstable ignition increases, leading to progressively more violent flashes up to a final 
flash with Helium shell mass w 0.02-0.1M Q . The high pressures of these last flashes allow the burning to 
produce the radioactive elements 48 Cr, 52 Fe, and 56 Ni that power a faint (M v = -15 to -18) and rapidly rising 
(few days) thermonuclear supernova. Current galactic AM CVn space densities imply one such explosion 
every 5,000-15,000 years in 1O U M of old stars (« 2-6% of the Type la rate in E/SO galaxies). These ".la" 
supernovae (one-tenth as bright for one-tenth the time as a Type la supernovae) are excellent targets for deep 
(e.g. V = 24) searches with nightly cadences, potentially yielding an all-sky rate of 1,000 per year. 
Subject headings: binaries: close — novae, cataclysmic variables — supernovae: general — white dwarfs 



1. INTRODUCTION 

Stellar evolution can lead to double white dwarf (WD) bi- 
naries, often in orbital periods (P or b) short enough that gravi- 
tational wave losses drive them into contact (Webbink 1984). 
Prior to making contact, these binaries are sources for space- 
based gravitational wave interferometers (Nelemans et al. 
2001, 2004) and have been found in numerous optical sur- 
veys (Marsh et al. 1995; see Nelemans et al. 2005 for the 
latest overview). For comparable mass WDs, the outcome 
at contact (when the lower mass WD fills its Roche lobe) is 
likely catastrophic (Webbink 1984; Marsh et al. 2004). 

For extreme mass ratios (hence the lower mass WD is pure 
He; Marsh et al. 2004), the mass transfer at initial Roche 
lobe filling is dynamically stable (see Deloye & Taam 2006 
for a recent discussion). Under the continued loss of orbital 
angular momentum from gravitational radiation (Faulkner et 
al. 1972; Tutukov & Yungelson 1979; Nelemans et al. 2001; 
Deloye, Bildsten & Nelemans 2005), the binary separates to 
P orb ps 60 minutes as the Helium donor gets whittled down to 
< O.O1M in over 10 9 years. This AM Canum Venaticorum 
(AM CVn) class of P orb < 60 minute binaries (Warner 1995; 
Nelemans 2005) has been known for 30 years, but recent har- 
vesting of the Sloan Digital Sky Survey (SDSS, Roelofs et al. 
2005; Anderson et al. 2005) and X-ray detections of possible 
direct impact binaries (Cropper et al. 1998; Israel et al. 2002; 
Ramsay et al. 2002; Marsh & Steeghs 2002) have greatly in- 
creased the sample. The number of AM CVn binaries found 
in SDSS imply a space density of (1 -3) x 10~ 6 pc" 3 , includ- 
ing those which started mass transfer 10 Gyrs ago (Roelofs et 
al. 2007b) 

We begin in $2] by discussing the fate of the Helium that 



accretes on the WD at rates 10~ 9 - 1O~ 6 M0 yr" 1 . The burn- 
ing is unstable for these accretion rates, leading to successive 
shell flashes likely to appear as classical novae. However, we 
find that the last unstable flash at low M's has a mass in the 
0.02-0.1M Q range, large enough to become dynamical and 
eject radioactive 48 Cr, 52 Fe, and 56 Ni. In §3, we show that 
the resulting thermonuclear supernovae are faint (My « —16) 
and rapidly evolving ".la" events. We close in §4 by deriving 
their occurrence rate in old stellar populations and discussing 
the prospects for their detection. 

2. HELIUM SHELL FLASHES IN AM CVN BINARIES 

Before discussing the Helium shell flashes expected from 
AM CVn binaries, it is important to mention the original work 
on thicker helium shell flashes in the context of Helium star 
donors (i.e. stars actively burning helium with initial masses 
M c / w O.4M ). These are studied as another potential pro- 
genitor for thermonuclear supernovae via the double detona- 
tion scenario (Nomoto 1982a, 1982b; Woosley et al. 1986; 
Woosley & Weaver 1994). Once the helium star comes into 
contact (Iben & Tutukov 1987, 1991), the mass transfer rate 
is set by gravitational wave losses (Savonije et al. 1986; 
Ergma & Fedorova 1990), leading to the M-Mj trajectory 
(from Yungelson 2007, this binary started at 80 minutes and 
came into initial contact after some He burning had occurred) 
shown by the dotted line in Figure 1 . These donors adjust their 
thermal state until P or b ~ 10 minutes (where Mj w 0.2M Q ), 
and then expand under further mass loss. Most of the Helium 
is added to the C/O WD at a rate of 2 - 6 x 10" 8 M Q yr -1 . 

The solid lines are Iben & Tutukov's (1989) estimated He 
ignition mass (Mj gn ) dependence on M for different C/O WD 
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FIG. 1. — Fully degenerate (semi-degenerate) Helium donor mass in AM 
CVn binaries are shown by the dashed line (long-dashed line from C. De- 
loye, priv. comm.) as a function of M, while the dotted line is for a He 
burning star (Yungelson 2007). The solid lines display Iben & Tutukov's 
(1989) estimates of Helium ignition masses for pure He accretion on WDs 
with M = 0.6,0.8, 1.0 & 1.2Mq as a function of M. The names of the two 
AM CVn binaries with M's in this range are shown at their respective M's. 

masses. Calculated with the He burning star scenario in mind, 
these are accurate for M < 1O~ 7 M yr -1 , and compare well 
to numerical work (Nomoto 1982a; Fujimoto & Sugimoto 
1982; Limongi & Tornambe 1991; Yoon & Langer 2004). 
This illustrates that the large He ignition masses (> O.2M ) 
needed to detonate the Helium and the underlying C/O are 
achieved in this scenario for O.6M WDs. Yoon & Langer 
(2004) found much lower Mi gn 's at these M's when they in- 
cluded large amounts of rotational dissipation at the depth of 
Helium burning. However, the work of Piro & Bildsten (2004) 
showed that the accreted material comes into co-rotation with 
the star at much lower depths, so that there is no rotational 
impact on Helium ignition. 

Our focus here is the Helium ignition expected in the AM 
CVn binaries that contain a non-burning helium donor and a 
temporally declining M. The short-dashed line in Figure 1 
shows the M-Mj relation (decreasing with time as the or- 
bit widens) for a fully degenerate He donor. We begin the 
plot for M < 2 x 1O~ 6 M yr" 1 , as the Helium burns stably 
at higher M's (e.g. for a > 0.2-0.27M Q He WD donor at 
P or b « 2.5-3.5 min; Tutukov & Yungelson 1996), potentially 
appearing as a supersoft source (van Teeseling et al. 1997). 
However, as M decreases, the Helium burns unstably via re- 
current flashes. For an accreting WD with M = 0.6M Q , there 
will be w 10 recurrent weak He flashes as the donor mass 
drops from 0.2 to O.1M . Deloye et al. (2005) and Deloye 
& Taam (2006) showed that a more realistic expectation is 
that the He donors have finite entropy (see also Roelofs et al. 
2007a), which increases (long-dashed line) the M for a fixed 
Ma, nearly to the dotted line case, which is the helium burning 
star donor after adiabatic expansion quenches the burning. 

In all these cases the "last flash" will have the largest M l%a . 
For simplicity, we use the intersection of M- lgn with the de- 
generate donor mass in Figure 1 to find the fiducial last 
flash masses M ign = 0.11, 0.09, 0.06, 0.04 & 0.02M Q for 
M = 0.6,0.8, 1.0, 1.2 & 1.36M WDs, shown as the short- 
dashed line in Figure 2. Given the uncertainties in binary 



10" 1 



s io- 



io- 



0.6 




1 

M core [M ] 



FIG. 2. — Minimum Mj gn needed to achieve ; nuc = fdyn and t„ 
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a function of the core mass (thick solid lines). An isothermal core of 3 X 10 7 
K is assumed (Bildsten et al. 2006). The dashed line shows the "last flash" 
Mj gn , with an uncertainty of 50% above and a factor of two below (hatched 
region). 

evolution, we allow for larger flash masses of 50% (the in- 
crease if we find intersections with the dotted line in Figure 
1) and lower flash masses of a factor of two (to accomodate 
the phase of the last flash and the possibility that M; gn may be 
lower than given by Iben & Tutukov (1989)), resulting in the 
hatched region of Figure 2. 

Once ignited, the flash proceeds by developing a convec- 
tive zone that moves outward as the base temperature rises 
due to the energy generation rate from the triple-a reac- 
tion, e nuc (Nomoto 1982a). The thermonuclear runaway can 
lead to a dynamical burning event and potential detonation 
(Taam 1980a, 1980b; Nomoto 1982b) if the thermonuclear 
timescale, f nuc = c P T / e nuc , becomes as short as the local dy- 
namical time, fdyn = H/c s , where H = P/pg is the pressure 
scale height and c s is the sound speed. 

In order to calculate these timescales, we integrate hy- 
drostatic equilibrium for an isothermal C/O core (including 
Coulomb corrections to the ion pressure) and a fully con- 
vective pure He envelope. Because M; gn can be a significant 
fraction of the core mass, the decreasing pressure during the 
convective evolution allows for hydrostatic expansion of the 
core, so we must calculate the structure of the core and en- 
velope consistently. (Our calculation thus differs from that of 
Fujimoto & Sugimoto (1982), who assumed that the core did 
not adjust during the flash.) For a given core mass and M; gn , 
the radial expansion of the burning envelope leads to a max- 
imum temperature that sets the minimum value of f nuc /fdyn- 
This minimum value decreases as M; gn is increased. Thus, for 
each core mass, there is a minimum M\ gn for a fully convec- 
tive, hydrostatic envelope that can become dynamical during 
its expansion. 

The thick solid lines in Figure 2 show the minimum Mj gn 's 
that achieve f nuc = fdyn and f nuc = 10fd yn (to allow for the pos- 
sibility of initiations in turbulent fluctuations; c.f. Wunsch 
& Woosley 2004). Also shown is the estimated range of 
the "last flash" M ign (hatched region). For a C/O WD with 
M > 0.7-0.9M Q , the last flash is likely to be dynamical. For 
lower mass WDs, whether or not the last flash is dynamical 
will depend on the binary evolution and accurate calculations 
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of M lga . However, stable mass transfer at Roche lobe con- 
tact favors extreme mass ratios, potentially leading to more 
massive WDs in AM CVn's (Marsh et al. 2004), as recently 
inferred (e.g. Roelofs et al. 2007a). 

These Helium flashes have ignition pressures lower 
(P pa 10 22 - 10 23 erg cm" 3 ) than previously studied in detail. 
Hashimoto et al. (1983) showed that constant pressure burn- 
ing at these P's can produce the radioactive elements 44 Ti, 
48 Cr, 52 Fe, and 56 Ni. We ran a few detonation calculations 
using a one-dimensional, explicit, Lagrangian, finite differ- 
ence scheme in which the hydrodynamics are coupled to a 
13 isotope a-chain reaction network. The network includes 
a-chain, heavy-ion, and (a,p)(p,j) reactions. Starting from 
the fully convective profile at the moment ? nuc = fd yn , we per- 
turb the temperature at the base of the envelope by 1%. A 
steady detonation forms on a dynamical time and propagates 
outward, steepening as it moves down the density gradient. A 
weak shock (over-pressure ~ 2) propagates downward into 
the C/O WD but does not ignite the Carbon at the enve- 
lope/core interface. In general, we find products of He en- 
velope burning similar to those found by Hashimoto et al. 
(1983). For example, for a 1 M Q C/O WD with a 0.06M Q 
Helium envelope, at the time the shock reaches the WD sur- 
face, 56 Ni is the most abundant ash (« 0.012M Q ), followed by 
52 Fe (0.0076M Q ), 48 Cr (0.007 1M Q ), ^Ti (0.005 3M ), and 
unburned helium (0.026M Q ). Further burning will occur as 
the matter cools, but we have yet to calculate it fully. 

We do not follow the shock deeper into the core of the WD, 
so the possibility remains for inwardly propagating shocks to 
detonate the Carbon at or near WD core (Livne 1990; Livne 
& Glasner 1991; Woosley & Weaver 1994; Livne & Arnett 
1995; Garcfa-Senz et al. 1999). However, our range of 
last flash M; gn 's is lower than every model that achieved a 
double detonation in the surveyed literature (Nomoto 1982b; 
Woosley et al. 1986; Livne 1990; Livne & Glasner 1990, 
1991; Woosley & Weaver 1994), save for one data point in 
Nomoto (1982a), which is at the upper bound of our range of 
last flash M ign 's. Note that Livne & Arnett (1995) achieved 
double detonations inside our range of last flash Mi gn 's, but 
they artificially induced He detonations in models that should 
not have become dynamical. Thus, it is likely that the last 
flash ignitions will not detonate the underlying WD, but this 
issue needs further examination. 

3. RADIOACTIVE OUTCOMES AND FAINT 
SUPERNOVAE 

Given the diversity of outcomes, we only present prelim- 
inary calculations of the peak luminosity and time to peak 
light, t„„ when M e = 0.02-0.1M Q of radioactive 56 Ni (decay 
time r = 8.77 d), 52 Fe (r = 0.5 d) or 48 Cr (r = 1 .3 d) is ejected 
from the WD surface at the velocity v w 15,000 km s -1 con- 
sistent with the nuclear energy release (minus the binding en- 
ergy as it leaves behind the C/O WD). The time of peak lu- 
minosity, r m , is when the age since explosion matches the dif- 
fusion time through the envelope (for an opacity k), yielding 
r 2 = nM e /7cv or 

( k \ 1/2 ( Me \ 1/2 ( 10 9 cm s-' \ 1/: 
r,„«2.5d^ 01cm2g _ I j ^ Q Q5M J ^ v J 

(1) 

where the factor 7 is found by both Arnett (1982) and Pinto 
& Eastman (2000). In Type la supernovae, the time to peak 
is 17-20 days (Conley et al. 2006; Garg et al. 2007), whereas 
the lower envelope masses in these explosions give r,„ ~ days. 
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FIG. 3. — The peak magnitudes versus time to peak light of faint AM 
CVn supernovae powered by 48 Cr/ 56 Ni (solid) or 52 Fe (dotted) for (left to 
right) k = 0.1,0.2, and 0.4 cm 2 g -1 . Each line shows the range of ejected 
masses from M e = 0.02-0.1M Q with velocities v = 12,000- 17,000 km s" 1 
(set by the underlying WD mass). The crosses show the 0.06Mq explosion 
for the same three values of k and for v ~ 15,700 km s -1 (brighter set) and 
v = 5,200 km s -1 (dimmer set). The long-dashed line denotes the region 
occupied by Type la supernovae, including sub-luminous 1991bg-like events. 

Though the rapid rise of these events presents an observa- 
tional challenge, our expectation is that they are very bright. 
This is because from Arnett's (1982) "rule" we know that 
the peak luminosity matches that produced instantaneously 
by radioactive heating. Hence the shorter r m allows us to 
see the power from the short-lived radioactivities, yielding 
peak bolometric luminosities of 1 - 5 x 10 42 erg s , some- 
what lower than Type la supernovae (Contardo et al. 2000). 
These events come close to matching the peak brightnesses 
of sub-luminous Type la supernovae like 1991bg, though the 
observed long decay time requires an ejected mass > 0.1M Q 
(Stritzinger et al. 2006). 

The radioactive heating from the 56 Ni decay chain is well 
known, but the radioactive heating from 52 Fe and 48 Cr needs 
amplification. The 52 Fe decays (via positron emission 55% 
of the time) with a half-life of 8.275 hours to the 1 + state of 
52 Mn, followed by a 168.7 keV 7-ray decay to the metastable 
2 + state 52m Mn, giving an average heating of 0.86 MeV. The 
52m Mn isotope has a 21.1 minute half-life and decays via 
positron emission to the 2 + state of 52 Cr, releasing 3.415 MeV. 
The 48 Cr isotope has a 21.56 hour half-life with nearly 100 
% electron captures to the 1 + excited state of 48 V, followed 
by a cascade that emits an energy of 0.42 MeV. The 48 V has 
a 15.973 day half-life, and an effective energy deposition of 
2.874 MeV. 

In order to estimate the appearance of these events, we cal- 
culated t,„ from equation (Q~|i with k's motivated by Type la 
SNe (Pinto & Eastman 2000; Mazzali & Podsiadlowski 2006) 
for the range of M e = M; gn of the dashed line in Figure 2. For 
the velocities, we assumed that the matter leaves the WD with 
a kinetic energy equal to that available once the nuclear en- 
ergy release unbinds the material. In order to calculate the 
absolute V magnitude, My, shown in Figure 3, we assumed 
the radioactive heating luminosity at t = r„, for ejected en- 
velopes of pure 56 Ni (solid curve), 52 Fe (dotted curve), or 48 Cr 
(solid curve). By coincidence, the radioactive heating curves 
of 56 Ni and 48 Cr are similar at 2-10 days. For the brightest 
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events ( 48 Cr and 56 Ni), 7k « 10,000-20,000K at maximum 
light (we assumed stellar colors to make our bolometric cor- 
rections), the ejected envelope is always radiation dominated, 
and the largest Mi gn 's are the brightest. The rapid decay of 
52 Fe means that the largest Mi gn 's are the faintest (and have 
?eff < 3,000K) since r m exceeds the 52 Fe half-life. The low 
Mgn 52 Fe events are brighter and have T e ff w 10,000 K. 

Clearly work remains, including complete ignition and det- 
onation calculations with full nucleosynthetic outcomes and 
resulting velocity profiles. In order to show what might be 
possible, the crosses detail the appearance of the O.O6M ex- 
plosion reported in §2 for the same three values of k, and two 
values of v. The lower value of v was considered just in case 
some of the nuclear energy in the helium detonation is lost to 
shocking the underlying C/O WD. Clearly, even in that case, 
these remain bright events worthy of observational interest. 

4. DISCUSSION AND OBSERVATIONAL 
CONSEQUENCES 

Our consideration of the maximum Helium ignition masses 
naturally attained in AM CVn binaries led us to detonations 
that primarily produce 48 Cr, 52 Fe, and 56 Ni. These rapidly 
decaying radioactive elements, combined with the low ejected 
masses, cause rapidly rising (2-10 day) events with My w— 16. 
Detection of such events would reveal the presence of AM 
CVn binaries in distant galaxies; all that remains is for us to 
estimate their rate. 

The measured local galactic density of AM CVn's 
(Roelofs et al. 2007b) yields a mass specific AM 
CVn birthrate in an old stellar population of (0.7 — 2) w 
10" 15 AM CVn yv- l LQ tK ~ l . If every AM CVn gives a visibly 



explosive last flash, this would give 1 event every 5000-15000 
years in a 10 n M Q E/S0 galaxy, 2-6% of the Type la SN rate 
in E/SO galaxies (Mannucci et al. 2005; Scannapieco & Bild- 
sten 2005; Sullivan et al. 2006). Since the time elapsed from 
reaching contact to the last flash is < 10 8 years, the occur- 
rence rate of these faint supernovae in an actively star form- 
ing galaxy depends on the rate at which double WDs come 
into contact (Nelemans et al. 2004). Additionally, if every 
AM CVn gives 20 weak Helium novae, then 1 in 5000-15000 
classical novae would be a Helium-rich event. 

For now, we use the local K-band luminosity density 
(Kochanek et al. 2001) to estimate the volume rate of 
.la's from the old stellar population to be (0.4 - 1) x 
10" 6 .Ia yr _1 Mpc -3 , 1 .6-4% of the local Type la SN rate. Up- 
coming optical transient surveys (e.g. Pan-STARRS-1, Pan- 
STARRS-4, Large Synoptic Survey Telescope) with rapid ca- 
dences should find these events. Observations with daily ca- 
dences to a depth of V = 24 would have an all-sky rate of up 
to 10 3 yr" 1 (typical distance would be « Gpc). Since these 
.la's evolve quite rapidly, frequent followups, such as those 
provided by the Las Cumbres Observatory Global Telescope 
(www.lcogt.net) will prove a valuable resource for this en- 
deavor. Their rapid evolution to small optical depths may also 
allow for strong 7-ray line emission. 

We thank D. Kasen and S. Kulkarni for discussions, J. Ste- 
infadt for rate calculations, C. Stubbs for naming these events 
.la's, and F. Timmes for providing the reaction network. This 
work was supported by the NSF under grants PHY05-51 164 
and AST02-05956. GN is supported by NWO Veni grant 
639.041.405. 
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